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Abstract 

We evaluate the 7r°7* — > 7 transition form factor and 77 decay widths for 
7r°,r/ c and rjb treated as qq bound states in the Bethe-Salpeter formalism, 
incorporating the dynamical chiral symmetry breaking and the Goldstone na- 
ture of the pion. In the chiral limit, the Abelian axial anomaly is incorporated 
analytically in our coupled Schwinger-Dyson and Bethe-Salpeter approach, 
which is also capable of quantitatively describing systems of heavy quarks, 
concretely rj c and possibly 7/&. 
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1. Recently, Bando et al. [|]] and Roberts have demonstrated how the Adler-Bell-Jackiw 
axial anomaly can be incorporated in the framework of Schwinger-Dyson (SD) and Bethe- 
Salpeter (BS) equations, reproducing (in the chiral limit) the famous anomaly result for 
7T° — > 77 analytically. This was also extended to the off-shell case 7*7r° — > 7, a process 
of renewed experimental interest because of the plans to measure it with high precision at 
CEBAF (see refs. in §,!). 

Refs. ||,[3| avoided solving the SD equation for the dressed quark propagator S by us- 
ing an Ansatz quark propagator. Then, thanks to working in the chiral limit, they also 
automatically obtained the solution of the BS equation: in this limit, when the chiral sym- 
metry is not broken explicitly but dynamically, and when pions must consequently appear 
as Goldstone bosons, the solution for the pion bound-state vertex rV, corresponding to the 
Goldstone pion, is given (see Eq. (Q), and, e.g., ||) by the dressed quark propagator S(q): 

S'\q)=A{ q 2 )i-B{q 2 ) . (1) 

Concretely, [fj,|3| used the solution for that is of zeroth order in the pion momentum p, 
and which fully saturates the Adler-Bell-Jackiw axial anomaly It is proportional to 

B(q 2 ) when (|I|) is the propagator of a massless quark, whereas its normalization is given |5j 
by the pion decay constant f n : 

T 7T (q;p 2 = M 2 = 0) = 2 l5 B(q 2 ) m=0 / f n , (2) 

(where the flavour structure of pions has been suppressed) . Thanks to using the chiral-limit 
solution (H|) (and to satisfying the electromagnetic Ward-Takahashi identity), |],|2| analyti- 
cally reproduced the famous axial-anomaly result independently of what precise Ansatz has 
been used [fj for S(q) and, consequently, independently of what interaction G^ u (k) has 
formed the 7r° bound state. 

In the present paper we argue that in the light of the recent experimental results on 
Vc — * 77 from CLEO [0], elements of these treatments concerning the electromagnetic 

interactions of dressed quarks, appear essential also for the understanding of the electro- 
magnetic processes of mesons in a totally different regime, far away from the chiral limit. 
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However, the above-described scheme which avoids solving BS equations PH§, is re- 
stricted to the chiral limit and very close to it. (Already when strange quarks are present, 
(0) can be regarded only as an "exploratory" Q expression.) This scheme obviously must be 
abandoned in the case of c- and 6-quarks, where the whole concept of the chiral limit is of 
course useless even qualitatively. Away from this limit, one must confront solving the perti- 
nent bound-state equation, which is determined by the interaction between quarks. On the 
other hand, one should also have the axial anomaly incorporated correctly. This consistency 
requirement is satisfied by the coupled SD-BS approach to qq bound states presented by 
Jain and Munczek PHIOH (and reviewed in a broader SD-BS context by |f^|II]), because this 
treatment in the chiral limit yields pions as Goldstone bosons of dynamical chiral symmetry 
breaking (D%SB). On the other hand, it has also successfully reproduced |8HT0|1 almost the 
whole spectrum of meson masses, including those in the heavy-quark regime, and also the 
leptonic decay constants (fp) of pseudoscalar mesons (P). We are therefore motivated to use 
this so far successful approach for calculating other quantities. In this paper we present the 
calculation of the vr°,^ c , and 7#, — > 77 decay widths, as well as of the 7*7r°-to-7 transition 
form factor. (Of course, in the case of 77^ we do not question the accuracy of the non- 
relativistic description, but want to see how well our approach does once the experimental 
data on rjb are obtained.) 



2. To define the interaction kernel for the SD and BS equations, [pHTOf use the modeled 



Landau-gauge gluon propagator given as the sum of the following two parts: 

CT(k) = [G uv (-k 2 ) + G m (-k 2 ))(g^ - ^) , (3) 

which define the separation of the propagator into: a) the well-known perturbative part, 
correctly reproducing the ultraviolet (UV) asymptotic behaviour unambiguously required 
by QCD in its high-energy, perturbative regime, and b) the nonperturbative part, which 
should describe the infrared (IR) behaviour. The infrared behaviour of QCD is however still 
not well understood, so it is only this latter, nonperturbative part that is in fact modeled. 
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From the renormalization group, in the spacelike region (Q 2 = —k 2 ), 

lQ7ia(0 2 ) —d ln[ln(x + j^-)] 

3 Q 2 g2i n ( Xo + ^_) i n ( Xo+ «_) 

where the two-loop asymptotic expression for a s (Q 2 ) is employed, and where d = 12/(33 — 
2JV», 6 = 2(3 2 /(3l = 2(19A 7 / /12 - 51/4)/(JV>/3 - 11/2) 2 , and JV> is the number of quark 
flavours. Following |ID[ , we set A 7 / = 5, Aqcd = 228 MeV, and = 10. For the modeled, 



IR part of the gluon propagator, we choose Gjr also from Ref. [JTO 

G IR (Q 2 )= 1 -^faQ 2 e-^\ (5) 

with their [JTOj parameters a = (0.387 GeV) -4 and fi = (0.510 GeV)" 2 adopted throughout. 

We obtain the dressed quark propagators ([1]) for various flavours by solving the SD equa- 
tion in the ladder (or "rainbow" ||TTJ) approximation (i.e., with bare quark-gluon vertices): 

s- 1 ( q ) = 4-m-zj 0- 4 YS(k)YG, u ( q - k) , (6) 

where m is the bare mass term of the pertinent quark flavour, breaking the chiral symmetry 
explicitly, (i.e., we do not use any Ansdtze for A(q 2 ) or B(q 2 ).) The case rh = corresponds 
to the chiral limit where the current quark mass m — 0, and where the constituent quark 
mass B(0)/A(0) stems exclusively from D^SB. For u and <i-quarks, assumed massless, solving 
of (|) leads to B(0)/A(0) = 375 MeV for the parameters from When m/0, the SD 
equation @ must be regularised by a UV cutoff A HJT0| , and m is in fact a cutoff-dependent 
quantity. In 0, A = 134 GeV, while m(A 2 ) is 680 MeV for c-quarks, and 3.3 GeV for 
6-quarks. This gives us the constituent mass B(0)/A(0) of 1.54 GeV for the c-quarks, and 
4.77 GeV for the 6-quarks. In the chiral limit, Eqs. (|l])-(|2D reflect the fact that solving 
of (|6|) for S(q) with m = is already sufficient to give us the Goldstone pion bound-state 
vertex (0) saturating the anomalous tt° — > 77 decay, and that is how our is obtained. Of 
course, we cannot avoid solving the BS equation like this for heavier qq pseudoscalars P, 
such as P = rjc, rjb- Their bound-state vertices Tp must be obtained by explicit solving of 

r P (q,p)=iJ 0^YS(q' + ^)T P (q\p)S(q' - ^)YG, u (q - q') , (7) 



the homogeneous BS equation again in the ("generalized" |2j) ladder approximation, con- 
sistently with (^). (For rj and rj', their mixing complicates the situation, so we do not solve 
for P — rj,r)' at this point.) Note that S is the quark propagator obtained by solving the SD 
equation (||) with the same gluon propagator G^. For pseudoscalar mesons, the complete 
decomposition of the BS bound state vertex T P in terms of the scalar functions Tf is: 

r P (q,p) = 75 { + j/Tf (q,p) + </T?(q,p) + & q] T^q,p) } . (8) 

(The flavour indices are again suppressed.) The BS equation (|7D leads to a coupled set of 
integral equations for the functions Tf (i = 0, 3), which is most easily solved numerically 
in the Euclidean space. Solving for T^ c , we also get the rj c mass M Vc = 2.875 GeV, whereas 
experimentally M Vc = 2.979 GeV. Nevertheless, this is just one example of the successful 
reproduction of heavy-meson masses in this approach |8|-|T0f. For r]},, where there are no 
experimental results yet, we predict M Vb = 9.463 GeV. 

3. The transition matrix element for the processes ir ,r] c ,rif, — > 77, and 7*7?° — > 7, is given 
by the famous triangle graph coupling two photons (with momenta k and k') to a neutral 
pseudoscalar P. Apart from the momentum-conserving delta-function ensuring p = k + k', 
and the photon polarization vectors e^(k, X),e u (k', A'), it is in essence equivalent to a tensor 
amplitude Tp U (k, k'), which is the time ordered product of two electromagnetic currents of 
quarks, evaluated between a pseudoscalar state \P) and the vacuum state (0|. By symmetry 
arguments, it can be written as pf T£ u (k,k') = e a(3 ^k a k l f5 T P {k 2 ,k' 2 ) , where T P (k 2 ,k' 2 ) is 
a scalar. When both of the photons are on their mass shells, k 2 = and k' 2 = 0, the decay 
width is T(P -> 77) = 7ra 2 Mp |T P (0,0)| 2 /4 . (P = tt°, rj c , rj b .) 

On the other hand, for the 7*71° — ► 7 transition, only the final photon 7 is necessarily 
on the mass shell, k 2 = 0, whereas k' 2 = —Q 2 7^ for the virtual photon 7*. The 7*7r° — > 7 
transition form factor F{Q 2 ) on the 7r° mass shell is then defined as 

P(Q 2 )=TX0,-Q 2 )/7>(0,0). (9) 
How to calculate Tp, or, equivalently, Tp V l The framework for incorporating electromag- 



netic interactions in the context of BS bound states of dressed quarks, advocated by (for 
example) [[l]-[§] and called the generalized impulse approximation (GIA) by is probably 
the simplest framework consistent with our coupled SD-BS approach. This means that in 
the triangle graph we use our dynamically dressed quark propagator S(q), Eq. ([[]), and the 
pseudoscalar BS bound-state vertex Tp(q,p) instead of the bare 75 vertex. Another essential 
element of the GIA is to use an appropriately dressed electromagnetic vertex q), which 

satisfies the vector Ward-Takahashi identity (WTI), 

(q'-q)^(q',q)=S- 1 (q')-S~ 1 (q) . (10) 

Our dressed quark propagator contains the momentum-dependent functions A(q 2 ) and 
B(q 2 ), so that the bare vertex 7^ violates ([H]), implying the nonconservation of the electro- 
magnetic vector current and of the electric charge. Therefore, as in GIA, we must also use 
such a dressed, WTI-preserving quark-photon vertex. Solving the pertinent SD equation for 
the true, dressed quark-photon vertex T M is a difficult problem which has only recently begun 
to be addressed |13[]. Therefore, it is customary to use realistic Ansdtze ||11|| . Motivated by 
its successes in we choose the Ball-Chiu (BC) jI3J vertex for Y^^q'^q): 

r»tf,q) = A + (q-,q^ + ^^{A„{q>\ q^-Q^ - B^q' 2 ,q 2 )} , (11) 

where H±(q' 2 ,q 2 ) = [H(q' 2 ) ± H(q 2 )], for H = A or B. Obviously, this Ansatz does not 
introduce any new parameters as it is completely determined for each flavour by the pertinent 
solution for the quark propagator ([!]). Its four chief virtues, however, are (i) that it satisfies 
the Ward-Takahashi identity (|10|), (ii) that it reduces to the bare vertex in the free-field 
limit as must be in perturbation theory, (Hi) that its transformation properties under Lorentz 
transformations and charge conjugation are the same as for the bare vertex, and (iv) it has 
no kinematic singularities. 

It is important to note that the correct axial-anomaly result cannot be obtained |2j 
analytically in the chiral limit unless a quark-photon-quark {qq^) vertex that satisfies the 
Ward-Takahashi identity is used (even if D%SB is employed and the pion does appear as a 
Goldstone boson, as comparison with |T5|JT6|] shows). 
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In the case of 7r° for example, the GIA yields the amplitude Tp U (k, k')\ 

T$(k, hf) = -N c J -J^tr{r^ - f , * + q - P -)S(k + q - f ) 

x r (k + q - \, q + V -)S(q + - |)} + (* <- tf, /* <- i/). (12) 

(The analogous expressions for r] c and r^, or any other neutral pseudoscalar, are straight- 
forward.) The number of colours N c arose from the trace over the colour indices. The 
u and d quark charges, Q u = 2/3 and Qd = —1/3, appeared from tracing the prod- 
uct of the the quark charge operators Q and T3/2, the flavour matrix appropriate for 7r°: 
tr(Q 2 r 3 /2) = (Q 2 U - Qj)/2. 

4. Let us first discuss the 7r°7* — > 7 transition form factor F(Q 2 ) for space-like transferred 
momenta. The existing CELLO data |T7[, displayed (following 0|) as Q 2 F(Q 2 ) in Fig. 1, are 
described well by a monopole curve with the slope at Q 2 = which is equivalent to the "in- 
teraction size" (r^o 7 )^ = 0.65 ±0.03 fm (defined in, e.g., [1 via (r 2 ^) = -6F'{Q 2 ) Q 2 =Q ). 
This was reproduced in various models. E.g., in the relativistic constituent quark model 



with an Ansatz for the pion wave function, Jaus [18 obtained a slope parameter equiv- 



alent to (r 2 ^^) 1 ^ 2 = 0.64 fm, although the transition amplitude at Q 2 = (ir° — ► 77 
transition), is 20% too small. In a generalization of the Nambu-Jona-Lasinio model, Ito 
et al. |I9| got F(Q 2 ) in excellent agreement with the data and with the monopole fit with 
(r^ 7r o 7 ) 1//2 = 0.65 fm. The amplitude obtained in the chiral perturbation theory gave (for 
small Q 2, s) F(Q 2 ) = 1 - (r 2 )Q 2 /6, with (r 2 „ ) 1/2 = 0.28 fm @], but introduction of vec- 



7"7T U 7/ 

')-'Vi /••(Y)- i = 1 

7?r 7/ I ' \ 77T 7/ 

tor mesons increased it to (r 2 ) 1 ^ 2 = 0.64 fm [21]. Bando and Harada's [22] Breit-Wigner 



forms for the propagators of p and u resonances, have infinite derivatives for Q 2 —>■ + , so 
that one cannot give (r 2 7r o 7 ), although their F(Q 2 ), when pure vector-meson dominance is 
assumed, almost coincides with the monopole curve with (r^ 07 ) 1/2 = 0.65 fm. (See Fig. 1.) 
We also display the Brodsky-Lepage interpolation curve |23| to the perturbative QCD fac- 



torization limit Q 2 F{Q 2 ) — > 8n 2 f n = 0.67 GeV 2 . However, our results (the solid line in Fig. 
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1) should in the first place be compared with F(Q 2 ) of the closely related approach of Frank 
et al. [0 (the dashed line). Our F(Q 2 ) is a little better, but it is in fact relatively similar to 
the curve of ||, considering that our B(q 2 ), obtained by solving of the SD equation (||) (in 
this case in the chiral limit), is very different from the Ansatz for B(q 2 ) used by ||. Also, 
our (r 2 ^^) 1 ^ 2 = 0.46 fm is practically the same as 0.47 fm of so that the "interaction 
size" discriminates even less between our respective S(g 2 )'s, which describe the structure 
of the pion @ in the chiral limit. It seems, therefore, that the differences in the internal 
structure between our pion and that of |§ play a relatively small role for 7r°7* — > 7, at least 
in the chiral limit. 

For the on-shell case ir° — > 77, the dependence on the pion structure falls out completely 
in this limit. Namely, our framework is in the chiral limit equivalent to [0-0, as demonstrated 
by the fact that, in this limit, with the solution (|j), we too can reproduce the famous 
anomaly result T.„-o(0, 0) = l/47r 2 f n analytically, in the closed form. The decay width is then 
given by r(7r° — > 77) = (a 2 /647r 3 ) (M^/f 2 ) flr2| , in excellent agreement with experiment 
- see Table I. Being a consequence of respecting the WTI ([IT]) and incorporating D^SB, 
this result is independent of our concrete choice of the interaction kernel and the resulting 
hadronic structure of 7r°, i.e., of B(q 2 ). This is as it should be, because the axial anomaly, 
which dominates tt° — > 77, is of course independent of the structure. It is then not surprising 
that calculations for 7r° — > 77 which rely on the details of the hadronic structure (instead 
on the axial anomaly) have problems to describe this decay accurately. 

Of course, the situation is very different for quark masses of the order of Aqcd and 
higher: only the numerical evaluation of Tp(0,0) is reliable in this regime. Moreover, the 
details of the chosen interaction kernel and the resulting propagator functions A(q 2 ) and 
B(q 2 ), as well as the bound state solutions, do matter in that regime (which gets further 
and further away from the domination of the axial anomaly with growing quark masses, as 
demonstrated vividly by the amplitude ratios in the last column of Table I). This is naturally 
the case with the heavy-quark composites i] c and 775. Their two-photon widths are also given 
in Table I. In the case of T]b, only predictions exist, as there are no experimental data yet. In 
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the case of r] c , our 77 width is (without any adjustment of model parameters) comfortably 
within the empirical error bars after the 1994 Particle Data Group average, and its 1995 
electronic update p4f , have been further updated with the new, 1995 CLEO |7|] result of 
(4.3 ± 1.0 ± 0.7 ± 1.4) keV. As opposed to this, it seems that most of the other theoretical 
approaches predict r] c — > 77 decay widths that will be either too large or too small after 
inclusion of the 1995 CLEO result in the Particle Data Group average. 

To elaborate on this, let us remark that the results for T(rj c — *> 77) in the constituent, 
nonrelativistic potential models, once in good agreement with experiment, have recently 
been shown pH] to rise to far too large values of 11.8 ± 0.8 ± 0.6 keV after the calculations 



have been improved by removing certain approximations. The estimates |25] of the relativis- 
ts corrections indicate that they can reduce the width back down to some 8.8 keV. This is 
not enough to bring the constituent quark model in agreement with experiment, although 
the relativistic corrections are substantial, corroborating the view that the relativistic ap- 
proaches to bound states, such as ours, retain their importance on the mass scale of the 



c-quarks. This view has recently been strongly advocated by a broad overview of the 
two-photon physics of mesonic gg-composites in the BS approach (but without D%SB and 
without a dressed qqj- vertex); it has concluded that the relativistic effects are important 
for any two-photon width, even for heavy quarkonia. It is interesting to note, however, that 
such BS calculations, typically yielding |26j T(r] c — > 77) below 4 keV, are consistent only 
with the 1995 CLEO result on T exp (r] c — > 77), whereas they are too low for all other |24j 



measurements, and even for the newly lowered average of 6.2 ±1.3 keV in Table I. {E.g., see 
p6| and other related BS calculations discussed therein. They underestimate 77-widths in 
the light sector even more.) 

The results of our approach are higher and in agreement with the empirical decay widths. 
Our results also reveal the importance of using appropriately dressed quark-photon vertices 
together with dressed quark propagators in yet another way, which has to do with internal 
consistency. To demonstrate the consequences of violating the WTI flH]|), and thereby of 
non-conservation of electric charge, we compare (in the fourth column of Table I) the decay 
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widths Y bare (P — > 77) obtained inconsistently, i.e., in the coupled SD-BS approach |L0| but 



using the bare electromagnetic vertex 7 M , with our results T(P — > 77) using the WTI- 
preserving dressed vertex (|TT|). While not so catastrophic as for the light 7r°, the violation 



of the WTI (|T0| ) is still very bad for rj c — > 77, and very noticeable even for r)b composed of 
very heavy 6-quarks, for which the dynamical dressing is not as important as in the light 
sector. 

The reader should note that we have not done any fine-tuning of the parameters, or 
of the gluon propagator form (§)-(|) which we used; these are the propagator and the 
parameters of Ref . |IJ , which achieved a broad fit to the meson spectrum and pseudoscalar 



decay constants. Therefore, in the consistently applied generalized impulse approximation, 
the coupled SD-BS approach genuinely, without fitting, leads to the adequate amplitude 
strength, which is otherwise too low. In other words, the BS approach j8|HT0f which is 
consistent with the ideas of |]-[3| on dressed quark-photon interactions, seems to be able to 
describe electromagnetic processes well even in the heavy-quark sector without fine-tuning 
of model parameters. On the other hand, since T(?7 C — > 77) does depend on the interaction 
kernel and the parameters determining the internal structure of i] c , making the measurements 
of the processes such as r] c — > 77 more precise can, through our theoretical approach, 
contribute to determining the nonperturbative gluon propagator more accurately. In the 
first place, this pertains to the infrared part, which is at present poorly known, but is also 
of significance for the determination of the QCD running coupling a s (contained in (|j), the 
ultraviolet part of our gluon propagator), for example, at the scale m c naturally sampled by 
the i] c system. 
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TABLES 



p r(p -> 77) pe^^p 77 ) 



^Wp _ 77 ) T^o(0,0) 

l(p E 77) r P (o,o) 



7T° 7.7 7.74±0.56 0.22 1 

t/c 5.3 x 10+ 3 (6.2 ± 1.3) x 10+ 3 0.43 3.95 

0.15 x 10+ 3 ? 0.73 133 



Vb 



TABLE I. Comparison of the calculated 77 decay widths (in eV) of TT°,r] c and r\h with their 
average experimental widths, where we have included the 1995 CLEO result @. The widths 
L feare (P — ► 77), obtained with the bare electromagnetic vertices 7^, are significantly reduced with 
respect to L(P — > 77) due to the non-conservation of charge even in the heavy-quark regime. At 
present, there are no experimental data on rjb, so the calculated mass of % had to be used in the 
phase-space factors, unlike 7r° and r\ c . 
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FIGURE CAPTIONS 



1: 7*7r — > 7 form factor. The solid line represents our results. The dashed line represents 



the results of Ref. [EJ. The dash-dotted line is the Brodsky-Lepage interpolation [23 



The line of open circles is the curve of |2^|, and open squares form the monopole 
curve corresponding to (r^o ) = 0.65 fm. Experimental points are the results of 
the CELLO collaboration |T7[|. 
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